Spin current in the Kondo lattice model 
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By using the projection operator technique it is observed that the strong Hund's rule coupling 
and s-d interaction in transition metal elements may lead to an effective coupling between the spin 
current and spin chirality. As a result, the spin chirality can be regarded as a driving force to 
produce a spin current. The spin current may give rise to a novel type of field acting on the spins. 
A spin battery is designed based on the interactions between the spin current and spin chirality. 
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Spin-dependent effects, such as giant or colossal mag- 
netoresistance and magnetization switching, arise from 
the interaction between spin of charge carrier and an ap- 
plied external magnetic field. In spintronics, instead of 
charge, electron spin carries information and it is possible 
that capability and performance can be enhanced in spin- 
tronics devices li 2 ^ Possible applications and fundamen- 
tal science involved make the study of spin-dependent 
transport an intensive field in condensed matter physics. 
Recently, several ways to produce a pure spin current are 
proposed, such as anomalous Hall effect in ferromagnetic 
metals 4 ' 5 , ferromagnetic resonance&i, and the spin cur- 
rent in a spin spiral stated* 9 - In this paper, we study spin 
current in the Kondo lattice model, a simplified model 
containing essential physics for many ferromagnetic met- 
als or semiconductors. We find a new type of interaction 
between spin current and spin chirality (defined later) in 
the strong Hund coupling limit. This indicates that the 
spin chirality can be regarded as a driving force to pro- 
duce spin current, similar as an electric field does to the 
electric current. Furthermore, the spin current induces 
a new type of an effective magnetic field. Based on the 
interaction between spin current and spin chirality, we 
propose a scheme to design a spin battery. As a demon- 
stration, we study the spin transport in the well known 
spiral state to bring out the above mentioned spin-current 
effects. 

The Hund's rule coupling and s-d exchange interac- 
tion are crucial in transition metal elements. It plays a 
key role in the formation of metallic ferromagnetism in 
colossal magnetoresistance materialsiSiii* 1 ^ and diluted 
magnetic semiconductors ^2*14 In this paper we focus on 
doped transition metal oxides and/or ferromagnetic semi- 
conductors. The minimal model for these materials is the 
Kondo lattice model, 



runs for the nearest neighbor pairs of lattice sites. In this 
paper we will set H = c — 1 and the lattice space is a and 
is set to unitii 5 " For the Hund's rule coupling, J is always 
positive. The model has been studied extensively^ The 
spin current operator along a axis in the tight binding 
approximation is defined as thl 



nn+a 



and 



nn+a 



2 y j {^n,a T <y<y' c n+a,a' c n+a,a T vv' c n,a 



(2) 

where a(— x,y,z) are vectors pointing to the near- 
est neighbor sites, which have one longitudinal and two 
transverse components. The spin polarized currents and 
electric current become equal in a fully polarized ferro- 
magnetic system. 

Before we discuss the spin transport of Eq.(l) in an 
arbitrary dimension, let us first to illustrate the physics 
by considering the well known spin spiral state in one- 
dimensional system. In general, the spin spiral state may 
not be the ground state of Eq.(l), however, it can be 
stabilized by external environments, such as by apply- 
ing external fields along different directions at the ends 
of the system^ Our purpose here is to demonstrate 
the consequence of the spin spiral state. We parame- 
terize S n by the polar angles 9 n and tp n , i.e., making 
the classical spin approximation. The spin spiral state 
is defined by S n /S = (sin 9 cos nip, sin 9 sin nip, cos 9) . 
The spin chirality defined as the cross times between the 
two neighboring spins: (S„ x S„+i) z = S 2 sin 2 9sinip. In 
the k- space, we introduce a spinor operator ^(k) = 

c la+<p/2,p c la-<p/2,l) • The mean ficld Hamiltonian is 
written as 



H = -t 



(n,m) ,( 



J 

2 



S„ 



(i) 



where c\ a and c n ^ are the creation and annihilation op- 
erators for conduction electrons, S n is the spin operators 
for the localized spin, and t are the Pauli matrices and to 
is the 2x2 identity matrix. The summation over (n, m) 



H = 1 £2&(k)H{k)$(k), (3) 



where H(k) = -g a T - g x r x - g z r z with g = u + 
2t cos ^ cos fca; g x = ^ sin 6; g z = -| cos 9 — 2t sin j sin ka. 
u is the chemical potential. The single particle Green 
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FIG. 1: Spin current via density of conduction electron in the 
spin spiral state of = it/2 and ip = 0.2n. 



FIG. 2: Spin current via spin chirality at J 
density of conduction electron 0.20. 



1.0 with the 



function in a 2 x 2 matrix form is given by 

G(k,iuj n ) = {iuj n TQ - ff(fc))" 1 

[(iu n + g )T - g x r x - g z r z 



(iuj n 



m) 



(4) 



where uj n = (2n + l)irkT with an integer n and T is the 
temperature. The two branches of spectra are 



oj± (fc) = —It cos ^ cos ka 



± 



■ COS( 



2i sin % sin ka) 



^ sin 2 i 



1/2 



(5) 



In the k-space, the expectation value of the longitudinal 
spin current operator is written in terms of the Green's 
function, 

< J *> = N i A + n F( UJ + - M) - A_n F {cj_ - (j,)] , (6) 



where 
A±(k) 



± sm — cos fca 
2 

i sin y sin 2 ka — ^ cos cos i| sin ka 
^ cos - 2t sin | sin ka) 2 + ^ sin 2 6* 

We also calculated the electric current I c z and found that 
I z — 0. Thus in the spiral state there exists a spin current 
without a charge current. In Fig. ^we plot spin current 
as a function of electron density in the state with 9 = tt/2, 
and ip = 0.27T for different J. We see that a larger J 
produces a higher spin current. For a finite J the current 
may change direction with change of density. In Fig. [5]wc 
plot the relation between spin current and spin chirality 
at the density of 0.2. The spin current increases with the 
spin chirality. 

In the large J limit, the problem can be simplified. 
An equivalent Hamiltonian of Eq.(l) or so-called double 
exchange model in the large J limit is given byi& 



E 

(n,m) 



(7) 



with t nm = t [cos cos ^Sr + sm ^ sm 



and a n = cos(0„/2)(l - n ni )c„ T + exp[iy>„] sin(0 n /2)(l - 
n ni)cni and satisfies the anticommutation relation, 



a k ,a' k , 



5k,k' ■ The energy excitations and their eigen- 



states can be obtained from the diagonalized Hamilto- 



-2t 



■ cos kc 



sm 



■ cos(fca + ip) 



a{a k , 



where a k is the Fourier transform of a n . In the case 
when the spins of conduction electrons are frozen by the 
localized spins, and the quasiparticles become spinlcss 
fermions. The two Fermi momenta are k±a = —k ± pn 
where fco = arctan[sin 2 9 sin if /(cos 2 9 + sin 2 9 cos ip)] and 
p is the ratio of the number of conduction electrons to 
the number of total lattice sites and p/a is the den- 
sity of conduction electrons, fco is nonzero when the or- 
der parameter for the spin spiralling is nonzero. Its 
lowest energy state is \GS) = f] a \ |0) ■ The 

spin current in the state can be calculated from I s z = 
J2 n (GS\ If m+1 \GS) /N. We have an analytical expres- 
sion for spin current, 



It = th 



sm pir 



'sin. ip 



1 -sin 2 0sin 2 ((/j/2) 



(9) 



The result is equal to what we obtain by taking the large 
limit in Eq.©. In the point of view of the Berry phase 
the spin spiral state acquires a constant Berry phase in 
the spin-dependent renormalized factor in the Hamilto- 
nian (Eq. [J, i.e., t nn+ i — |t nn +i | e - ^. The non-zero 
phase 5tf drives a spin current. 

By analogy with electric field and resistance, we intro- 
duce the concept of spin current resistance or its recipro- 
cal, spin conductance. The role of S„ x S m here is analo- 
gous to the local electric field. In the spin spiral state the 
spin chirality is non-zero, (S„ x S n +i), = S 2 sin 2 9 sirup. 
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In Eq.@ the current is proportional to the spin chiral- 
ity (S n x S n _)_i) , and can be written in a compact form 
l s z = c ■ (S„ x S n+ i) z /S 2 by introducing a spin conduc- 
tance 



sin pir I 



>sin 2 (>/2) 



(10) 



in a large J limit. The spin conductance is determined 
by the filling of electrons and the magnetic structure. 
Just like an electric conductance, the spin conductance 
reflects intrinsic properties of spin transport in materi- 
als. It depends on the spin-dependent scattering mech- 
anism of electrons. In Eq. i|10[l . for the case of p — > 
or — > 1, c is proportional to the density of electrons p 
or the density of holes 1 — p. For a small spin chirality 
sin 2 9sm 2 (tp/2) — > 0, the spin conductance aprroaches to 
[sinpTrJ/Tr, which is proportional to the averaging kinet- 
ical energy of the charge carriers in this tight binding 
theory. From the numerical results in Figs. 1 and 2, we 
deduce that the spin conductance also depends on J. It 
becomes more complicated for the finite J cases. Simi- 
larly, we can also introduce a spin voltage, which is also 
proportional to the length along the spin spiraling path, 



V, 



spin 



= J (S„ x S n +i) ■ dl/S 2 where dl = r n+1 - r„ 



IS 



the difference of two position vectors of neighboring spins. 
V sp i„ is an important parameter characterizing the abil- 
ity to produce a spin current. 

Thus far we have introduced several important con- 
cepts, such as spin chirality or spin current resistance, 
through the example of the classical spiral state. To de- 
termine the microscopic mechanism to produce a spin 
current, we go back to the Kondo lattice model (Eq.(l)) 
in the large Hund's rule coupling limit and address the 
quantum spin case. The following discussions are not 
limited to one-dimensional systems and the conclusion is 
independent of dimensionality According to the Hund's 
rule coupling, the energy of empty or double occupancy 
of conduction electrons is zero. For a single occupancy 
the conduction electron and localized spin can form a 
spin 5 + 1/2 state with — J5/2 and a spin 5 — 1/2 state 
with +J(S + l)/2. In the large J limit the spin 5-1/2 
state and double occupancy should be excluded as they 
have much higher energy than the spin 5+1/2 state. 
This process can be realized with the help of the pro- 
jection technique^ For that, we introduce the spinor 

. The dressed oper- 
are used here to 



ators (1 



r )cjj a , instead of c\ 



avoid the double occupancy due to the strong coupling. 
The projection operator for the spin 5 + 1/2 state is de- 
fined as 



K = <t>iK , 



(ii) 



where P+ = [S„ • r + (5 + l)r ]/(25 + 1). The operator 
■p+ forces the electron spin at site n to be fully polarized 
with the localized spin at the same site. The projection 
operator for non-occupancy is V!^ = (1 — n„^)(l — n n> i). 



The total projection operator to exclude the spin S- 
1/2 single occupancy and double occupancy is V = 
EI C^n + T^n ) • Therefore the equivalent Hamiltonian in 



the strong Hund's rule coupling limit is expressed a; 
H eff = VHV = -t£ [0+P+P+> m + h.c] , 



,12.20 



(12) 



where a constant is omitted in the last step. In the large 
S approximation the Hamiltonian of fEa ll2fl becomes the 
double exchange model in Eq.(J7J. Using a mathematical 
identity 

(S„ • r) (S m • r) = S„ • S TO r + i (S„ x S m ) • r, (13) 
the Hamiltonian can be rewritten as 21 

(5 + i) 2 i 



H = 



(25+1) 
t 



]T1 X! fc " 



(25 4 


I) 2 


(5 + 


l)t 


(25 4 


I) 2 







~~ "2 ^ ^ (Sn S m ) ' ]VI nm 

nm 

J2(SnXS m )-K m , (14) 



(25 +1) 2 

where k nm = <p+T cj) m + <p+ n T cf) n and M nm = c/)+T(j) m + 
0+r0„. The first term is the kinetic part for conduc- 
tion electrons excluding double occupancy, like the t- 
tcrm in the t — J model. The second term contains 
an effective Heisenberg-type spin-spin exchange interac- 
tion. Approximately, the effective coupling is Jferro = 
—t {k nm ) /(25 + l) 2 , proportional to the kinetic energy 
in the sense of the mean field approximation. The kinetic 
energy —t (k nm ) in the tight binding approximation is al- 
ways negative if the electron filling p < 1 such that the ef- 
fective coupling is ferromagnetic. This is consistent with 
the double exchange picturei^SiSi The third term is the 
interaction between local spin S„ and the vector M„ m . 
For a paramagnetic state, (M nm ) is equal to zero, but 
in a fully polarized state, its value is equal to (k nm ) , i.e. 
proportional to the kinetic energy. Among the four terms 
the most important observation is from the last term: it 
contains an interaction between the spin current operator 
If lm and the spin chirality S„ x S m . By analogy with the 
coupling between the electric current and electric field it 
is straightforward to understand why a non-zero spin chi- 
rality induces a pure spin current. Although we cannot 
derive an analytic expression for a finite J, it is expected 
that the mechanism to produce a spin current remains 
the same. 

A new effective field and spin battery. To clarify the 
physical meaning of the last term, we come to calculate 
the rate of change of a local spin. From classical mechan- 
ics or from the Heisenberg equation of motion the rate of 
change of angular momentum S„ is equal to the torque 
p n x which acts on the spin, 



dS n /dt = -gp B S r . 



x K ff , 



(15) 
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Interface 

FM Metal Spin Spiralling Material 




FIG. 3: Schematic view of DC spin battery. Arrows in the FM 
metal as a spin source indicate the magnetization is perperdic- 
ular to the interface. The spin spiraling material provides a 
non-zero spin spirality to drive spin current along the axis. 



where the effective field experienced by the spin is 

[knn+sSn+8 + {S + l)M„„ + j + I„ n+S X S n+ $] . 

The first two terms, again are related to particle's orbital 
motions. The third term is the torque between the spin 
current and the local spin. A new type of force or ef- 
fective field is caused due to the interaction between the 
spin current and local spin or local magnetization. The 
force may be responsible for magnetization precession by 
the spin injection into a magnetic conductor— and the 
anomalous Hall effect in ferromagnetic materials Fur- 



ther discussion along this line will be published elsewhere. 
Our observation will lead to potential applications. For 
example, by means of the spin spiraling properties we 
can design a spin battery to produce a pure spin current. 
The battery contains two parts as shown in FigEI The 
main part is a spin spiraling material. The spins spirals 
about the axis which is perpendicular to the interface 
with a ferromagnetic metal film such that the spin volt- 
age along the axis is non-zero. The spin chirality in the 
material will produce the spin accumulation at two termi- 
nals. The magnetization of the film is also perpendicular 
to the interface. The film is a source of spins which will 
flow out of the battery driven by the spin chirality. The 
key element in this battery is the spin spiraling materials. 
Possible realization of spiral states have been discussed 
by several authorS ) 27 ' 28 i 29 i 30 and the state may also be 
stabilized with help of the external magnetic field and 
exchange coupling between localized spins. 

In conclusion, we have shown that the strong coupling 
between the spins of charge carriers and local spins leads 
to an effective interaction between spin current and spin 
chirality in the system. Non-zero spin chirality may give 
rise to a pure spin current and in return a spin current 
can cause a novel type of force on the spins. We confirm 
these findings by studying in detail the spin transport in 
the spin spiral state. 

The work is supported by grants from the Re- 
search Grants Council of Hong Kong, China (Project 
HKU7088/01P) and US-DOE. 
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